Calculation of the interspecies s-wave scattering length in an ultracold Na-Rb vapor 
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We report the calculation of the interspecies scattering length for the sodium-rubidium (Na-Rb) 
system. We present improved hybrid potentials for the singlet X 1 E + and triplet a 3 E + ground 
states of the NaRb molecule, and calculate the singlet and triplet scattering lengths a„ and a t for 
the isotopomers 23 Na 87 Rb and 23 Na 85 Rb. Using these values, we assess the prospects for producing 
a stable two-species Bose-Einstein condensate in the Na-Rb system. 
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I. INTRODUCTION 

The s-wave scattering length a plays a central role in 
the description of atom-atom collisions at ultralow tem- 
peratures (T <C 1 mK). In this regime, the cross section 
for elastic collisions, <7 e ; ~ 7ra 2 , and the cross section for 
inelastic spin-exchange collisions, <r ex ~ ir(ctt ~ a s) 2 are 
both expressed in terms of a ,1]. The scattering length 
is also a critically important parameter in the physics of 
Bose-Einstein condensates (BECs) . For a bosonic atomic 
species i, a BEC is stable only if <Zj > 0. In addition, ef- 
ficient evaporative cooling demands that a e i 3> o~ ex 0- 
A priori calculations of a are thus of fundamental in- 
terest, and the quest for BEC in alkali-metal atoms has 
spurred on efforts to calculate the scattering length in 
many atomic species. 

A majority of the work on scattering lengths has con- 
centrated on interactions between like alkali- metal atoms. 
However, the recent production of Bose-Fermi mixtures 
in 6 Li- 23 Na 0, 6 Li- 7 Li HQ, and 40 K- 87 Rb || and 
dual-species BEC in 41 K- 87 Rb has renewed interest 
in binary mixtures of ultracold gases. In these systems, 
the interspecies scattering length a\2 is the basic quan- 
tity parametrizing the interactions between component 
atomic species, ayi determines the efficiency of sympa- 
thetic cooling en route to the formation of quantum de- 
generate mixtures, and in the case of two-species BECs 
(TBECs), a\2 determines the stability and miscibility of 
these mixtures @- However, relatively few calculations of 
interspecies scattering lengths exist (3, This is due, 
in part, to the incomplete characterization of diatomic in- 
teraction potentials for many pairs of alkali- metal atoms. 

In this paper we calculate the scattering lengths for the 
Na-Rb system. We construct the NaRb potential from 
a combination of spectroscopic data and precise long- 
range interaction parameters, and use a simple method 
for calculating the singlet and triplet scattering lengths 
for the isotopomers 23 Na 87 Rb and 23 Na 85 Rb. The Na- 
Rb system is interesting in part because its an obvious 
candidate for TBEC. Both of the component species have 
been condensed and the condensates have been studied 
in detail The Na-Rb TBEC has been treated theo- 
retically, and its properties are sensitive to the value of 
the interspecies scattering length a^a-Rb 0- 



This paper is divided into four sections. In Sec. [HI 
we discuss methods for calculating aNa-m, and address 
the source of errors involved in these calculations. In 
Sec. IIIII we introduce two hybrid potentials for the NaRb 
molecule, and contrast our potentials with previous re- 
sults. In Sec. II VI we present our calculations for a,Na,-Rb, 
and discuss the feasibility of producing a TBEC in the 
Na-Rb system. 

II. SCATTERING LENGTH 

If the potential V(r) is known for all r, then the scat- 
tering length a can be calculated 0. The procedure 
is to numerically integrate the radial Schrodinger equa- 
tion for low collision energies to large values of r. At 
large r the numerical wave function u(r) is "matched" to 
an asymptotic form x(r) by requiring that ^lnu(r) = 
J^\nx(r) at the match point r — r m . The asymptotic 
wave function can be written as 

Xi=o(r) ~ r [sm(Ar) + tan(<5 ) cos(fcr)], (1) 
k 

where k is the atomic wave vector and 5 is the s-wave 
phase shift. In the limit that k — ► 0, the wave function 
asymptotically approaches a straight line as r — > 00 . The 
scattering length is given by the r intercept of this line 
|13| . and can be formally defined in terms of the s-wave 
phase shift as 

a = _ Urn M*0 . (2 ) 

For low collision energies, r m should be large to ensure 
that the numerical wave function attains its asymptotic 
behavior. The total integration time can be reduced by 
calculating the corrections to the asymptotic wave func- 
tion at smaller values of the match point r = r m . For 
example, Marinescu has shown 0] that the wave func- 
tion for large r may be written as x(r) = cte a (r)-\- j3e.p{r), 
where the functions e Q , !l a(r) are solutions to the differen- 
tial equations e^^(r) = [2fj,V(r)/h ]e aj p(r) subject to 
the boundary conditions e a (r) — ► r and ep(r) — ► 1 as 
r — > 00. These equations cannot be solved exactly. How- 
ever, if an analytic expression for the long-range potential 
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is known, then the functions e a .p(r) may be estimated to 
arbitrary precision using a method of successive approx- 
imations. The scattering length is given by a — —[3/ a, 
and can be found by applying the usual continuity con- 
dition at r = r m . 

The uncertainty in the calculated value of a depends 
upon the reliability of the potential V(r). Gribakin and 
Flambaum |la | have shown that for diatomic potentials 
which fall off at long range as —C^/r 6 , the scattering 
length is given by 

a = a o [l-taa{$(J5 = 0)-7r/8}], (3) 

where a a is a "mean scattering length", and Q(E) is the 
semiclassical phase, defined as 

where ri nner (r outer ) is the inner (outer) classical turning 
point of the potential at energy E and \i is the reduced 
mass of the colliding atoms. As can be seen from Eq. 
the scattering length is infinite if = 0) = (n — 3/8)tt, 
where (n=l,2,3,...). This situation occurs if the nth vi- 
brational state of the potential V(r) is barely bound at 
E = 0. In general, V(r) will not admit a barely bound 
state. However, variations of the potential within its es- 
timated uncertainties will shift the energies of its bound 
states, and states lying closest to dissociation experience 
the largest shifts. A bound state may even be intro- 
duced or removed from the well, depending on the size of 
the potential shift and the proximity of a bound or vir- 
tual level to the dissociation energy. As states are added 
or removed from the well, the scattering length passes 
through ±oo. Therefore, if the interatomic potential is 
not known well enough to predict whether or not a barely 
bound state exists, then a cannot be specified within fi- 
nite bounds. 

Because of the extreme sensitivity of a to the bind- 
ing energy of the highest vibrational state of the inter- 
atomic potential V(r), the most precise calculations of a 
in alkali-metal atoms typically rely on the spectroscopy of 
bound states near dissociation. Two-color photoassocia- 
tion or Raman spectroscopy is used to resolve these lines 
to high precision 0] . In the absence of near-dissociation 
spectroscopy, the scattering length may still be calcu- 
lated, but the accuracy and precision of such a calcula- 
tion is limited by the accuracy and precision of the inter- 
atomic potential V(r). In the case of alkali-metal dimers, 
spectroscopy is sparse and near-dissociation spectroscopy 
is nonexistent. However, the potential V(r) may still be 
"assembled" from RKR (Rydberg-Klein-Rees) [lTj data 
and well-known analytic expressions for short- and long- 
range potentials to create a "hybrid" potential valid for 
all r. In assigning error bars to a, care must be taken 
to ensure that variations of the hybrid potential within 
its estimated uncertainties do not introduce or remove 
bound states from the well. For example, uncertainties 
in the Cq coefficient of NaK allowed for additional bound 



states in its hybrid potential, and frustrated attempts to 
determine the scattering length for some isotopomers . 
Fortunately, this is not the case for NaRb, as shown be- 
low. 



III. NaRb POTENTIALS 

Compared to many alkali- metal dimers, the ground 
states of the NaRb molecule are relatively well-known. 
Rovibrational states to within 5% of dissociation have 
been observed in both the triplet a 3 E + and singlet 
X 1 T, + potential wells [l^. |20| . The rotationless inter- 
atomic potentials V s (r) and V t (r), corresponding to the 
X 1 Yi + and a 3 S + states, respectively, have been deter- 
mined through RKR analysis, and a direct fit to the sin- 
glet spectrum [2(3 using a modified Lennard-Jones (MLJ) 
|2l| parametrization has also been performed. Ab initio 
ground state potential curves for NaRb have been calcu- 
lated j^. However, these curves are not very accurate, 
so we do not use them in constructing our interatomic 
potentials. 

At large r, the NaRb interatomic potential is accu- 
rately represented by a sum of two independent contribu- 
tions, the exchange and dispersion energies. The disper- 
sion energy is given by a well-known expansion in powers 
of r _1 : 

V ^p( r ) = -[^6+^ + ^w)- ( 5 ) 

The coefficients C n may be calculated from a knowledge 
of atomic polarizabilities |23| . The exchange interaction 
is calculated u sing the surface integral method of Smirnov 
and Chibisov |2J|, which yields 

V ex {r) = ±^J{A,B,a,(3,r) r ^+i-^m- 1 e -(<^P)r ) 

(6) 

where a 2 /2 and f3 2 /2 are the ionization energies (in 
atomic units) of each atom, and r is assumed to be in 
units of Bohr radii. The function J(A, B, a, (3, r) can be 
expanded in a power series ^2 n (Jnr n (a — (3) n )/n\ whose 
coefficients J„ are expressed as integrals that must be 
solved numerically. The complete long-range potential is 
then given by 

V LR (r) = V disp (r)±V ex (r). (7) 

The exchange energy is positive (negative) for the 
triplet (singlet) state. As r — * oo, the long-range in- 
teraction potential is dominated by the well-known van 
der Waals potential —Ce/r 6 . The exchange interaction 
is expected to become important inside the LeRoy ra- 
dius RheRoy |25j |. beyond which the potential is well ap- 
proximated by the dispersion energy alone. For NaRb, 

RheRoy ~ 11 A. 

The NaRb molecular potentials can be modelled by 
smoothly joining RKR data to the long-range interaction 
potentials. Zemke and Stwalley (Z-S) have constructed 
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TABLE I: Our chosen values for the parameters of the NaRb 
long-range potential Vlr{t). 



TABLE II: Scattering lengths found from our hybrid poten- 
tials for 23 Na 85 Rb and 23 Na S7 Rb, in units of Angstroms. 



Parameter 



Value 



C 6 a 
Cs b 
Cio b 

OtNa 

Pm 

ANa d 
Rb 

Jo e 

Jl E 
J 2 C 



1.293xl0 7 
3.4839 x10 s 
1.1552x10"' 
0.61459 
0.55409 
0.76752 
0.56945 
1.4197xl0~ 2 
6.0963xl0~ 4 
1.9537xl0" 3 



a C% given in units of cm -1 A . See Ref. 1271 , 

b Cg and Cio g iven in units of cm -1 A 8 and cm -1 A 10 , respec- 
tively. See Ref. pi . 

c The quantities are expressed in atomic units. See Ref. |29j. 

''The constants A and B are related to the size of the wavefunction 
of each atom in the region of interaction. See Ref. \'M\ . 

e We found that the exchange energy was adequately represented 
in our region of interest by the first three terms of the expan- 
sion J(A, B, a, /3, r) = '^2i n (JnT n {a — /3) n )/n\. Here we use atomic 
units. See Ref. 151. 



such hybrid potentials for the a 3 E + and AT 1 £ + states of 
the NaRb molecule |2a. More complete spectra of the 
NaRb singlet state ;20) and a more precise estimate of the 
Cq coefficient for NaRb have since become available, 
allowing us to construct new hybrid potentials for both 
the NaRb X 1 E+ and a 3 E + states. For r < 11 A, we use 
the recent MLJ potential to model the state. Our 

potential for the a 3 £ + state is identical to the Z-S triplet 
potential for r < 13.5788 A. Our long-range potential 
differs from that used by Z-S in a number of ways. We 
use the Cs and Cio dispersion coefficients recommended 
by Marinescu and Sadeghpour 28], but choose for Cq 
the highly precise value calculated by Derevianko et al. 
[23 ■ For the exchange energy, we used the heteronuclear 
expression of Smirnov and Chibisov given by Eq. 
Both the singlet MLJ and triplet RKR potential curves 
are joined smoothly to our long-range potential as given 
by Eq. (JJJ. Our complete long-range potential is given 
in Table |U 



IV. RESULTS AND DISCUSSION 
A. Na-Rb scattering lengths 

To calculate the singlet or triplet scattering lengths, we 
choose the hybrid potential V s (r) or Vj(r) and integrate 
the radial Schrodinger equation at E — from r inn er to 
the match point r m using the Numerov algorithm [3_l|. 
The reduced mass \x is given by M^ a M.Bb / {M^ a + Mm) , 
where M labels the atomic mass of either 23 Na, 85 Rb, or 



Isotopomer 



as 



ObT 



23 Na 85 Rb 
23 Na 87 Rb 



167+ 50 

10 ' -30 

551? 



59t 9 2 



51 



+9 



87 Rb. We expect the isotopic correction to the internu- 
clear potential V(r) to be negligible [2J}. Following Ref. 



[14|. the scattering length is given by 



ue. 



u e c 



ep - u'ep 



(8) 



where u(r) is the numerically integrated wave function, 
and the functions e a ^(r) are determined from the long- 
range potential, as discussed in Sec. [H] The primes de- 
note derivatives with respect to r. We found that a 
fourth-order approximation to e Q (r) and ep(r) guaran- 
teed convergence to a reliable value of a at a match point 
r m = 100 A. Our results are summarized in Table IfTI 

Because the scattering length is very sensitive to the 
details of the interatomic potential V(r), it is important 
to ensure that our calculated values of a are stable with 
respect to changes in V Sj t(r) within their known exper- 
imental or theoretical uncertainties. These include un- 
certainties in the value of the dissociation energy D e , the 
inner and outer turning points of the RKR potentials, the 
binding energy of the observed vibrational states, and the 
coefficients of the long-range parameters. We estimate 
our errors by calculating the change in the semiclassical 
phase $(£' = 0) due to the error in each parameter of our 
potential. These "phase errors" are then summed in 
quadrature to give a total phase error A4> to tai- We con- 
vert this value into a scattering length error using Eq. 
©■ 

For the triplet state, the error in a was estimated 
with respect to the uncertainties in D e , C 6 , and the 
RKR turning points. In the case of the singlet state, 
the MLJ parametrization allowed us to simultaneously 
vary all parameters in a statistically meaningful way. We 
treated the MLJ parameters as random variables with a 
well-defined mean and standard deviation. The phase 
$>(E = 0) was then calculated for 100 "random" poten- 
tials. The phase error was determined by examining the 
distribution of phases. This phase error was used to cal- 
culate the error in a s . 

As shown in Tabled our error bars are small. This is, 
in some ways, a fortuitous result. Had a bound or virtual 
state been closer to dissociation, variations in the poten- 
tial may have caused the phase to pass through a region 
where a — > oo. Because the scattering lengths are rela- 
tively small, they are more stable with respect to changes 
in the corresponding potential. In addition, our hybrid 
potentials are reasonably well constrained. This is due, 
in part, to the observation of bound states relatively close 
to dissociation in the NaRb triplet well, which enabled 




TABLE III: Recently determined scattering lengths for 23 Na, 
85 Rb, and 87 Rb, in units of Bohr radii. 



9 10 
r (Angstroms) 

FIG. 1: A plot of the difference between the Z-S and MLJ 
X E + potentials. 



Z-S to reduce the uncertainty in the dissociation energy 
to A(D e ) = ±0.1 cm -1 , and the very precise calculation 
of the Cq coefficient by Derevianko [23 ■ 

One interesting consequence of the new singlet poten- 
tial is the appearance of an additional bound state. We 
found that our complete potential for the X 1 !] 4 " state 
supported 83 bound states, whereas the corresponding 
Z-S potential only supported 82 bound states. This can 
be understood by examining Fig.^ which shows the en- 
ergy difference between the Z-S singlet and MLJ poten- 
tials for 6 A < r < 11 A. Because the MLJ potential 
is everywhere deeper, the wave vector k(E = 0) will be 
larger for all r. In semiclassical terms, the wave func- 
tion will build up more phase $ in this potential. In our 
case, $mlj — §z-s ~ so that another bound state 
appears in the MLJ molecular well. The energy differ- 
ences between the two curves is due primarily to errors in 
the extrapolation procedure used by Z-S to connect the 
short- and long-range components of their hybrid singlet 
state potential. 



B. Two-species Na-Rb BEC 

Having calculated the Na-Rb scattering lengths, we 
now consider the properties of a mixture of Na and 
Rb condensates. The interactions in a Na-Rb TBEC 
are parametrized by three scattering lengths: a^ ai aub, 
and aNa-Rb- The values of a R b and a^a are known 
to high precision, and are given in Table IIIII We con- 
sider a Na-Rb mixture in which the Na atoms are in 
the state \Fpf a ,mN a ), and the Rb atoms are in the 
state |-FRb, rriRb}, where F is the total angular momen- 
tum of the atom and m is its projection on the quan- 
tization axis. We represent the two-atom state by the 
kct \Fjsr a ,TriN a ; FRb,mnb). To calculate the scatter- 
ing leng ths we use the low-energy elastic approximation 
|34l l35j , which assumes that elastic collisions dominate 
the total cross section for collisions. This approximation 



Species 


a s 


at 


23 Na a 


19.1 ±2.1 


65.3 ±0.9 


85 Rb b 


970.5+420 


-388 ± 3 


87 Rb b 


90.4 ±0.2 


98.98 ± 0.04 



"See Ref. Bl. 
b See Ref. |3- 

gives for the scattering lengths 

a = a s P s ± a t P t , 



(9) 



where P s and Pt are the probabilities of the atoms be- 
ing in a singlet or triplet state, respectively. To calcu- 
late the probabilities P s and Pt, we project the state 
\FNa,m Na ; F R b,m R b) onto the states \S,m s ;I,mi), 
where S and / refer to the total electronic and nuclear 
spin of the two-atom system, respectively, while ms.i are 
their projections onto the quantization axis. This basis 
is useful for characterizing the system at smaller internu- 
clear distances where the exchange energy dominates. In 
this region, F^ a and F R b are no longer "good" quantum 
numbers, and the singlet and triplet states are labeled by 
S = 0, 1, respectively. 

To calculate the projections C^'™^f Rb ' mRb = 

(S,ms;I,mi\FNa,mN a ;F Rb ,mRb), w e perform the an- 
gular momentum recoupling of the four quantum num- 
bers Sno,, In a, Sub and I R b by making use of the 
Wigner 9-j symbols and standard Clebsch-Gordan al- 
gebra. We calculated the complete recoupling matri- 
ces Usi-ff for both 23 Na 85 Rb and 23 Na 87 Rb, and ex- 
tracted the probabilities P s and Pt for all input channels 
\FNai m Na'-i FRbi m Rb) ■ The scattering length for an ar- 
bitrary input channel is then given byEq. ©. 

Knowledge of a enables us to calculate the cross section 
for elastic collisions, a e \ = Air a 2 . Elastic collisions me- 
diate the rethermalization of atoms during evaporative 
cooling and sympathetic cooling. Furthermore, knowl- 
edge of a s and a t allows us to characterize the inelastic 
losses in the system. The dominant two-body mechanism 
for the loss of atoms from a trap is spin-exchange colli- 
sions. In this type of collision, the internal spin states of 
one or both of the atoms changes. In the case of magnet- 
ically trapped samples, such collisions can cause atoms 
to be ejected from the mixture by sending them into non- 
trappable spin states. More generally, the atoms may be 
ejected if the spin reorientation energy is converted to 
kinetic energies greater than the depth of the (magnetic 
or optical) potential confining the mixture. In the elas- 
tic approximation we can write the cross section for such 
inelastic processes as [35| 



Mifir{at - a s y 



(10) 



where Mif is a factor that depends on the asymptotic 
hyperfine states involved in the collision. Letting primes 
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denote the asymptotic output channel, we have 

-1 2 



M if = 



X! (Cs=oC's=o - C s =iC' s=1 ) 



(ii) 



where C is the projection coefficient defined above, and 
the indices have been suppressed. 

The achievement of a miscible two-component BEC 
places a number of constraints on the three relevant 
scattering lengths. Efficient sympathetic cooling re- 
quires a large magnitude of a^a-Rb- Collisional sta- 
bility against spin-exchange collisions requires small val- 
ues of <r ea; , which implies that the difference between a s 
and at be small. Dynamical stability of the individual 
BECs requires ajva > and am > 0. In the Thomas- 
Fermi approximation, the criteria for stability implies 
the existence of a critical value of |a./v a -_Rb| = a c above 
which the two-species condensate cannot coexist. The 
criteri a is given by \a Na _ Rb \ < a c = j^a Na a Rb , where 
7 = VM N ~M Rb /(M Na + M Rb ) ||. For a Na ^ Rb < -a c 
the attraction between the condensates overwhelms the 
repulsive interaction within each condensate and they 
collapse, while for aNa—Rb > &c the mutual repulsion 
of the two condensates is too great for them to overlap 
at all. 

The single-species scattering lengths are positive for 
both 23 Na- 23 Na and 87 Rb- 87 Rb collisions, which allows 
for single-species BEC in either atomic species. Using 
Eq. JHJ and the scattering lengths given in Tables HT1 and 
IIIII we calculate a c for all asymptotic states in a 23 Na- 
87 Rb mixture. The near equality of a s and at implies 
that <i e \ will be approximately the same for all states, 
and that <j ex will be small. Therefore interspecies elastic 
collisions will dominate inelastic spin-exchange collisions. 
Using the mean values for the triplet and singlet scatter- 
ing lengths given in Table [n] we find that for all asymp- 
totic two-atom states, a > a c . Taking into account the 
known uncertainties in the various scattering lengths, we 
find that the inequality a > a c still holds. We therefore 
conclude that a stable, miscible TBEC in a 23 Na- 87 Rb 
mixture is not possible. 



Next, we consider the mixture 23 Na- 85 Rb. This is an 
interesting case, since BEC has only been observed in 
85 Rb by utilizing a Feshbach resonance to tune the scat- 
tering length of the \F = 2, m F = -2) state To sim- 
plify our analysis, we eliminate from consideration those 
states for which the 85 Rb single-species scattering length 
is negative. Of the remaining states, we choose states 
that are lossless with respect to both homo- and het- 
eronuclear spin-exchange collisions. Because of the large 
positive singlet scattering length in 85 Rb, there is a large 
variation in a e i from state to state. If we use the mean 
values for the interspecies scattering lengths given in Ta- 
ble^ we again find no asymptotic states that satisfy the 
condition for TBEC stability. 

V. CONCLUSIONS 

We have derived hybrid potentials for the X 1 E + and 
a 3 £ + states of the NaRb molecule. We compare them 
to other recently derived potentials, and we discuss why 
our potentials are preferred. We have calculated the sin- 
glet and triplet scattering lengths from these potentials 
for both 23 Na 85 Rb and 23 Na 87 Rb. Using the elastic ap- 
proximation, we have calculated the scattering length for 
all two-atom asymptotic hyperfine states for both iso- 
topomers. The cross sections for elastic and inelastic 
spin-exchange collisions can be found using these values. 
Applying the Thomas-Fermi approximation criterion for 
TBEC stability, we find no two-atom asymptotic states 
for which a NaRb TBEC is stable. Further experimen- 
tal studies of ultracold Na-Rb vapors, including efforts 
to produce a TBEC in the NaRb system, will help re- 
fine our knowledge of the interatomic potentials and test 
these conclusions. 
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